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SUMMARY 



In the present report a.n: investigation is made on a 
flat plate in a two-dimensional compressible flow of the 
effect of compressibility and heati-ng on the turbu-lent 
frictional drag coefficient in the boundary layer of an 
airfoil or wing radiator. The analysi-s is based on the 
Prandt 1-Karm^n theory of the turbulent bo.undary layer and 
the St odola-Cr occ 0, theorem on the linear- relation betv/een 
the total energy of the flow and its velocity. (See ref- 
erences 1 and 2, ) Formulas are obtained for the velocity 
distribution and ;the frictional dra-g law in a turbulent 
boundary layer with the compressibility effect and heat 
transfer taken in-to ac-count.- 'rt, i:s found that with in- 
crease of compressibility and temperature at full retarda- 
tion of the flow (the temperature when the velocity of the 
flow at a' g"iv-eri poirit is r educ ed' "to ■ 2iero- in case of an 
adiabatic process in_ the gas) at a constant ' R^j. , the 
frictional drag" "c o'ef f ic'i e'nt c^ de"'cre'a'=!es ,.. both of these 
factors acting in the same sense. 

.INTRODUCTION -- 



In the present .pape.r an attempt -is made t.o generalize 
the theory of the Prandtl-Karman turbulent boundary layer 
to the case of flow of a compressible gas at high veloci- 
ties in the preg-ence of a temperature gra'dient from the 
wall to the gas. Sin-ce the theory itself ("mixing path" 
theory) is partly -.em.pi-ric-al , the generalization suggested, 
whi ch' *pTOb-a-bly c o-n-t--ai-n.3- ^an-. a-r.b4 t.rary,.-.-e;le.men.t,,- regui^res 
experimental conf i-rmat i on Moreover some error is. 
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introduced "by c e rt ai n , a s s'amptri on s that ar-e made to simplify 
the problem. The principal of th-ese assumptions are the 
following: ' 

(1) In the fundamental differential eouation for the 
velocity, di st rihiit i on in the TD.ound ary ."laye.r. : 



the frictional shear stress t, having a dependence on 
y that is not initially knbwn, if? to a first ap.proxiraa- 
tion eauated to its value at the wall, although with in- 
creasing distance from the wall it actually tends to zero. 
To increa'se the -a'ccuracy it is posit^ihle to proceed as 
follows: Havi^hg found'the Tplociiy d i p t rihut i on in the 
boundary layer by the method given, t is determined from 
the dynamic eaui li'brium taking account of the force's of 
inertia. From-the preceding eouation the velocity distri- 
bution is d.etermined to a" second ppproicima.tion and -so on, 
the process, how'ev-er, being, very complicated. 

(3) In th-fe" lamin.ar boundary layer the' Prandtl number 

Pr = ' i,Sj ,BqvLat,e.pl to u'nit.y although for pListom.ic gases 

it is actually eaua'l to 0.8. The error' arising from this 
assumption is s'ni'ali", "however. ' ' - • ' ■ 

(3) The nondimensional thickness of the laminar sub- 
layer and also the jump in the velocity gradient in pass- 
ing from laminar to turbulent, which, according to the 
Karman theory, are absolute constants, vjeve likewise re- 
garded as constant although they may vary in the present 
Case. The variation cannot, however, be theoretically 
determined. ' • ' 

(4) The condition- of-'the .transition from the pure 
laminar to the turbuTent --lawyer is arbi t r ary^- and is intro- 
duced only to Identify' the pr^blein' in question. Actually 
this transition depends on t.h.e xurbulence of the basic? 
flow. On account of th-i's indet ermrnat.enes s a Gerta.in ad--^ 
ditive constant .enters in. the formula for Re as a func- 
tion of Of o.f the order of ,10^ . (according to tests on 

liauids). . ' . . - 
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The theory proposed 'gives somewhaiJ more accurate re- 
sults than the Karman theory in which the resistance law 
ohtained for small v'eloclties is' maintained the same, suh- 
stituting only for p. and p their values at the wall 
hut not taking account of the effect of the compressihility 
on the velocity profile in the "boundary layer. The final 
answer can only he given, of course, "by experiment. 'Qual- 
itatively "both theories arrive at the same result r- namely, 
a decrease in Cf with increasing compressihility or heat- 
ing of the wall at constant Re. 



I, rUNDAMEHTAL D I FFEEENTI AL EQ,UATI01I 
■ INTSaEATION - BOUNDARY CONDITIONS 



Assume a thin smooth r e ct angiil ar plate situated in a 
plane parallel stream of a compressihle gas flowing with 
large constant velocity past a thin, smooth rectangular 
plate set at zero angle of attack. At the same time the 
plate is warmed by heat supplied from some heat source 
and' "by the friction so that a difference in temperature 
arises hetween the plate and the flow and heat is trans- 
ferred from the plate to the flow. 

The investigation is directly concerned with the 
prohlera of the effect of large velocities and of heat 
transfer on the frictional drag in the turbulent bound- 
ary layer about a plate on the supposition that this 
layer is very thin and has a quite definite boundary. 

Assume that th.e x-axis in the flow direction is 
along the plate, and the y-axia in the perpendicular di- 
rection from the , boundary layer toward the free flow. 
The variable velocity within thp boundary layer is denoted 
by u. Moreover, with the usual notation: 

p density 

viscosity 

T f rioti.o^n.al ..shear stress 

T absolute -temp'eirature 

0 specific heat at constant pressure 
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c^' specific heat kt constant .volume 

Cp/°y = k ■ • 

The c OTresponding varia"bles in the free flow' will "be ' 
denoted by '• 

p, U- , T; at the wall hy p;^, ui*. T*. 

The pressure p, as usual in the theory of the boundary 
layer, is assumed as constant. 

The n ondi men s i onal velocitipp and the magnitudes 
characterizing the state of the gas are functions of non- 
dimensional coordin5ites and therefore depend on certain 
parameters which are n ondi men s i o.ial c on:h inat i ons of scale 
ratios, end conditions of the plate, and physical con- 
stants of the gas; the shear stress t and the unit heat 
flow q and also the magnitudes p, p*, T^^ (.pressure p 
assumed constant) giving the state of the gas at the 
plate serving as end conditions. The density and temper- 
ature and T^ ar'= taken as scale ratios for p and 
T. The scale of velocity is expressed by the ratio 



T 



and the length scale as the ratio 



~ = — . ■ (2) 

where p.* is the value of the viscosity coefficient M> 
at the plate. 

With the introduction of these scales the parameters 
of the solution can only appear as nondimensional combina- 
tions of physical constants and end conditions. 

In the absence of compressibility and heat . transfer 
such combinations do not exist, as a result of which 
there is obtained, as' is 'known, a universal velocity dis- 
tribution not depending on any parameters, For this 
reason these scales are hereinafter -regarded as "univer- 
sal." 
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Then 



n = — n . ^^. 1 ,.. (3) 



\ 

u 

cp = — 



the square of the Bairstow num'ber: 



(k - 1) JcpT^ 



where U is the velecity in the free Gtream, 



Moreover 



V °f* 



wher e 



T ! 



and 



/ 



(4) 



a = -S:- = if^- (5) 



s 

U 3 , , 

= B (6) 



k = r:- = 1.4 (7) 



a « ii_--J^ b'^, (8) 
2 ^ 



(9) 
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1 ~ TT- ~ «- = "J (10) 

is the nondi mensi onal temperature of the flow hrought to 
rest adia"batically . 

It is assumed that in the e ompres sihle gas the futida- 
c^ntal differential equation of the turhulence theory of 
irarman remains valid; namely (reference 3), 




where _ is an ahgolute constant equal to 0„4 and t 
the frictional shear at a disXiance y from the plate. 

'fthen the coordinate sr (the distance from the lead- 
ing edge of the plate j is sufficiently large, however, 
the inertia forces may he neglected. Because of this and 
from the constancy of the pressure p it follows that T 
changes slowly v/ith y^ In ^_^irs,t_ approximation its 
value at the wall is assumed for 

Expressed in the universal s\3ales, equation (11) 
as Slimes the form 

In order to ohtain the differential eauation for cp 

it is only necessary to express o" terms of cp . 'This 

can he done with uhe aid of the known fact of the linear 

dependence of the total energy <-^c^T + u^/2 on the ve- 
locity u, (See lefeiencos 1 and 2,) 

In aondj mens i onal repre s en :ati on : 

1 = 1 ~ fr, ^4 + a, (13) 
h^ h / 
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or 



1 — 



(13*) 



Substituting equation (13*) In equation (IS), extracting 
the square root of both sides of the eauatlon, and using 
the minus sign in front of the radical, sinoe the expres- 
sion on the left in negatire, leaves 



V 



1 + 



4a 



(14) 



Setting 



(v < 1} and 



2 



V 



1 + 



4a 



(15) 



= A. 



(16) 



affords 



(17) 



whence, after integration 



ln-^ = X cos"^ v-{-lnC 



or 



a?] 



(18) 



I 

I 
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At the meeting of the laminar su'blayer and the turbulent 
layer 

where 

5 = 11.5. ^ (20) 

The determination of the houndary conditions proceeds 
from the asaumption of linear distribution of u in the 
laminar sublayer 

( reference 4-) 



or 



and hence 



? = >i (21') 
c?(s) = s. (22) 



i'urther, setting 



r V --^1 



results in 



fVa 



(23) 



, — ( cos"' Vt — COS"^ f ) 

(24) 



A second integration with the given end conditions 
gives 

where 



^The magnitudes s and f to a first approximation are 
assumed constant, their valuee being taken from Karman. 



HACA Technical Menorandum No. 1053 



9 



Or, if the original magnitudes are used again: 

1 



Tj=S-|- 




X 



— COS 



*^ h ^2 



(27) 



Expanding the expression on the right in and 



2 •/~ci 



and retaining terms through the second order approximates 
to 



, / , , lU a 3 10 2 \ 



I -^^^^ I Xtf -I-' 

^ 2z y ^ 

-■-ij, (28) 



where 



X'f = W 
XS = K'l 



(29) 



Ae a result of the substitution of numerical values, equa- 
tion (28) hecomes i icooo\i 

ttS 1.15+ -^(-g g — hzf- 16.222 ) + 



•0,034662 e'^ 



(30) 



As a— >0 and <0~>0 in equation (28), there is obtained 
in the limit 

hence 



9 = s + ^ln ((r,-s)x/+ll 



(31) 



or, let 
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-L- s = 5+-!^^ = ^ (32) 

then e<|uatlon (31) nay "be written in the forn 

,p = i-ln(7i + a) + &, (31') 

thus the well-known logarithmic law is ohtained. 
Further, for 

J, = 3, - (33) 

where 6 is the boundary layer thickness, or in nondimen- 
sional form, 

= = (35) 

From equation (28) is found the relation 

«./¥-+^{[>-T-i+.(F+T)_|/i + 

II. DETBHMINATION OF THE DRAG AT EACH 
POINT OF THE PLATE 

The inte-gral relation of Karman reads: 

s 

d r ..,r, (1) 
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Let 



(2) 



On introduction of the universal scales formula (l) 
aeeumee the form 



dRx . = A= 



(3) 



The integration of equation is carried out for 

the laminar sulilayer and the turbulent layer (fig. 4) (on 
the assumption that the turbulent layer begins at the 
point where the thickness of the laminar sublayer is 
equal to that of the purely laminar layer). In the lam- 
inar sublayer 



that is, 



'f = 



(4) 



s 

I' 



T( (A — -ri) df] 



1 — 



a-rj- 



A J 




' \ ^ h j 



, (5) 



where the first integral corresponds to the purely laminar 
layer, the second to the laminar sublayer of the turbulent 
layer, and the third to the turbulent layer. The first 
two integrals are small, as will be shown. 

tid in: 



It resultl 



d9~~ f 



■t.h 



^ h 2 T^a 



1 + 



Aa 



(6) 



\ J- 
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Developing the right-hand side of equation (6) in a 



(0 



Taylor series in </a and to terms of the second 

order approximate^ to 



1 *(-f-^) 

=7' 



d7i 1 * - 1 ' + 6l2" + + 4A 

1 + 



^ (»-^ + S'f + S=) + ^^('^ + s) 



6A' 



(7) 



or, from the original notation 



or 



« «' — ifi . 

ATI e ' , , a w „ 



62 



(70 



Substituting this expression in the last integral of 
equation (5) and carrying out the integration gives 



~ X'' I / tjZ \ 3 ^ 2 



z' + [ 1-f a 



CO j 

) + «(«',='-3)} + 



^-2 + 



(8) 



where L Is the value of the last two Integrals in equa- 
tion (5) for 2 3 wj (the lower limit) and 



e'dz 



(9) 



HACA ^•chn^cal Meaorandua Ho. 1053 



13 



1« the well-known integral exponential function for which 
tallies are available (for example, Jahnke and Snde). 'or 
■t equal to-16, 17, and 18 the valueu of. It^(t) were 

computed with the aid of the expansion of the function 
fiz) a in a McLaurin series which was then integrated 

term "hy term: 

" ln2r + z + -;5^ + ^^+ ... (10) 



/ 



2-2! 



3-3! 



By^ talcing 15 terms in equation (10), an accuracy to the 
fifth decimal place is assured. Terms of the form f"/n ni 
were computed with the aid of 10-place logarithm tables 
(which assure the required accuracy). The following re- 
sults were obtained: 

z=15 £',(15) = 234 956 
2=16 £(16) = 595 560 
z = 17 £-,(17) = 1 516 637 
2 = 18 £",(18) = 3 877 898 

By substitution of numerical values, formula (8) be- 
comes 



■a — cu 



■j 4 — |- a + 2.79io |z-f-|G—15.556a+ 18.160 10} + -^-^!"- 



— 2.708 (14.223a ^- 10.58(.j) £",• (z) + 39.0625 



1181.901a 



+ (7.748a — 19.403(1)) (2 In z — 1) — (15.693 + 2.1 69u)) z -\- 
+ (^ + ^ 0) + ~ aj. 447.746 + 41.981 +231. 866a + 42.498Trj 



(8') 



If the small terms not having e^ and S^^ ( z ) as factors 
are disregarded, it affords the approximation 
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S 1 .354 e^{z' - 42+ 6) + 1.354 a + _ 
_f'2 ^a_2.79u)^z — (15.556a— 18.160to) + 



3 ^ ; ^ ^3 

+ 85^336a j _ ^ .^^^ (i4.223a + 10.58«)') {z). (1 1) 

xUp* 

From the values of R_» ■ computed by the preceding 

xUp 

formula the values of H^^ = referred to the magnitudes 

u 

in the free flow were computed from the relation 



or, since 



and 



hence 



With the aid of the relation 



1 — ^ a = uj 

■' * 



there is ohtained 



Rx = 



[1 - (a + <i^)3'^2 
and - 

log H,^ « Ig Rjc* - - log [1 - (a + to)] (12) 

2 
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^ c 



f T' 



[1 - (w + a)] (13) 



The values of logS^ and were computed for: (l) 

a, = 0, W = 0; (8)03 = 0, a= 0,05, and a = 0 . 1 ; (s) 
Ct = 0,'co = 0.Q5, and 0) = 0.1; (4) a = 0.05, CO = 0.05, 
and to = 0.1; and (5) a = 0. 1, 00 - 0. 05, and 0) = 0.1 
in the interval from z = 7 to z = 18. Por these val- 
ues of a and (o the velocity distributions weye found 
for 2=7, 12, and 18, These results are given in 
tables 1 and 2, and shown graphically in figures 1 to 6 
(notation according to equation (39)). From inspection 
of figures 3 to 6 it is seen that, with increasitig 
and CO for the same » coefficient Cf decreases 
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TABLE I.- VELOCITY DISTRIBUTION 



w 


u 


H for 


01=0, 


»=0.0S 


« for 


<» = 0; 


i>=0,l 


u for 


a = 0, >i>=0.05 


u for 


a = 0; 


".=0,1 


u for a 


=0,05; M ~ 0,05 


a for a 


= 0.5; 


«> = 0,1 


ufor^ 


= 0,1; u 


= 0,05 


II for 0 


= 0,1; 


» = 0,1 




<u = 0 


z = 7 


z = 12 


2=18 


2 = 1 


j= 12 


z= 18 


2 = 7 


z= 12 


z = 18 


2 = 1 


z= 12 


z= 18 


2 = 1 


1 

2=\2 


z= 18 


z = l 


z= 12 


2=K 


2 = 1 


z= 12 


z = 18 


2 = 1 


1 = 12 


2= 18 


0 

b 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 


6,29 
15,13 
39,16 
104,45 
282,0 
764,6 
2076 
5642 
15336 
41686 
113312 
30S011 


6,27 
15,24 
40,07 
109,2 
303,2 
850,7 
2407 
6862 
19693 
56876 
165197 
482284 


6,28 
15,17 
39,47 
106,1 
289,2 
794.0 
2189 
6057 
16819 
46854 
130966 
367304 


6,29 
15,15 
. 39.30 
105,2 
285,2 
777,6 
2127 
5827 
15995 
43982 
121164 
334377 


6,26 
15,36 
40.98 
114.0 
324,3 
936,7 
2738 
8081 
24049 
72065 
217093 
656525 


6.28 
15,20 
39.78 
107,7 
2U6.4 
8-'3,2 
2301 
6472 
18302 
52023 
148619 
4265S6 


6,29 
15.16 
39,43 
105,9 
288,5 
790,7 
2177 

eon 

16653 
4C279 
129005 
360712 


6,25 
15,25 
-40,24 
109,8 
304,(5 
847,6 
2,172 
6658 
18727 
52765 
148891 
420682 


6,27 
15,20 
39.79 
107,5 
294,9 
813.1 
2249 
6235 
173i4 
48147 
134048 
373741 


6,28 
15,18 
39,58 
106,5. 
290,6 
796,9 
2192 
6037 
16S68 
45992 
127147 
351811 


6,21 
15,37 
41,31 
,115,1 
326,0 
930,6 
2668 
7674 
22119 
6,3845 
18447! 
533352 


6,25 
15,27 
40,41 
110,7 
307,7 
861,5 
2422 
6828 
19292 
54612 
154818 
439470 


6,26 
15,22 
40.00 
108,0 
299,2 
823,2 
2307 
6433 
18400 
50302 
140982 
395640 


11,37 
23,34 
79,16 
217,8 
606,2 
1697 
4778 
13519 
3841? 
109640 
314087 
902965 


11.40 
29,22 
78.11 
212.6 
583,1 
1607 
4437 
12291 
34132 
95000 
'265013 
741044 


11,42 
29,18 
77,73 
210,7 
574,8 
1573 
4318 
11863 
32863 
89973 
248310 
686187 


11,33 
29,46 
80,23 
223.2 
628,1 
1780 
5074 
14535 
41811 
120720 
349667 
1015635 


11,38 
29,29 
78,73 
215,7 
595,8 
1654 
4610 
12884 
36110 
101465 
285783 
806773 


11,40 
29,27 
78,15 
212,7 
583,3 
1606 
4433 
12259 
34394 
94283 
262145' 
730016 


11,36 
29.46 
«0,07 
222.6 
627,2 
1783 
5109 
14738 
42775 
124829 
36,5983 
1077206 


11,40 
29,25 
78,42 
214,1 
590,2 
1633 
.4549 
12706 
35615 
100169 
282666 
800337 


11,42 
29,19 
77,86 
211,3 
578,0 
1586 
4368 
12048 
33520 
92270 
256151 
712522 


11,32 
29,58 
81,14 
228,0 
649,1 
1866 
5405 
15754 
46167 
135909 
■401563 
1189876 


11, 3J 
29,32 
79,04 
217.3 
603.0 
1684 
4722 
13299 
37593 
106634 
303436 
866065 


11,40 
29,23 
78,28 
213,4 
586,6 
1619 
4483 
12443 
35052 
96580 
269986 
756351 



TABLE II.- VALUES OF LOG Rx AND Cf REFERRED TO THEIR FREE STREAM VALUES 



2 


1 = 0 


<u = 0 


0 = 0, i» = 0,05 


ci = 0, w = 0,l 


a = 0,05; i» = 0 


3 =0,1, 111 = 0 


d = 0,05 


111=0,05 


a = 0.05, 111 = 0,1 


a = 0,1, 


«) = 0.05 


a = 0,1, 


i« = 0,1 


log/?^ 








log/?. 








log/?. 




log/?x 




log/?. 




log/?. 




log/?. 




7 


4,60304 


0,00653 


4,661 


0,00620 


4.720 


0,00588 


4,656 


0,00620 


4,711 


0,00588 


4,715 


■ 0,00588 


4.77J 


0,00555 


4,771 


0,00555 


4,832 


0,00522 


8 


5,18576 


0,00500 


5.250 


O'.00bl75 


5.315 


0,00450 


5,253 


0,00475 


5,302 


0J)0450 


5,309 


0,00450 


5,372 


0,00425 


5,368 


0.00425 


5,433, 


0,00400 


9 


5,74784 


0,00.i90 


5.818 . 


0,00376 


5,888 


0,00356 


5,810 


0,00376 


5,872 . 


0,00356 


5,880 


0,00356 


5',950 


0,00337 ' 


'5.942 


0.00337 


6,012 


0,00317 


10 


6,29411 


0,00320 


6.370 " 


0,00304 


6,445 


0,00288 


6,360 


0,00304 


6,425 


0,00288 


6,135 


0,00288 


6,510 


0,00272 


6,500 


0,00272 


6,574 


0,00256 


11 


6,82794 


0.00264 


6,910 


0,00251 


6,989 


0,00238 


6,898 


0,00251 


6,965 


0,00238 


6.976 


'0,00238 


7,054 


0.00224 • 


7,044 


0,00224 


7,123. 


0,00211 


12 


7,35175 


0,00222- 


7439'- 


0,00211 


7.522 


0,002(10 


7,424 


0,00211 


7,496 " 


0.00200 


7,510 


■0,00200 


7,592 


0,00189 


7,580 


0,00189 


7.661 


0,00178 


13 


7,86735 


0.00189 


7.960 


0.00180 


8,048 


O.OOl^O 


7,943 


0,00180 


8.019 


O.0O17O 


8,033 


0,00170 


8,119 


0,00161 


8,106 


0,00161 


8.191 


0,00151 


14 


8,37603 


0,00163 


8.473 


0,00155 


' 8,565 


0,00147 


8,456 


0,00155 


8,534 


0,00147 


8,551 


0,00147 


8,638 


0,00139 


8,623 


0,00139 


8.712 


0,00130 


15 


8,87903 


0,00142 


8,981 


0,00135 


9,076 


0,00128 


8,962 


0.001J5 


9,043 


0,00128 


9,060 


0,00128 


9,152 


0.00121 


9,136 


0,00121 


9,226 


0,00114 


16 


9,37699 


0,00125 


9,48 1 


0,001 19 


9.584 


0.00112 


9,464 


0,00119 


9,546 


0,00112 


9,565 


0,00112 


9,660 


0.00106 


9,643 


0,00106 


9,736 


0,00100 


17 


9,87065 


0,00111 


9,983 


0,00105 


10,085 


0,00100 


9,')60 


0,00105 


10,046 


0.00100 


10,065 


0,00100 


10,163 


0,00094 


10,148 


0.00094 


10,242 


0,00089 


18 


10.36053 


o.oooyi 


10.477 


0.00094 


10.584 


0.00089 


10.454 


0.00034 


10,542 


0,00089 


10,564 


0.00089 


10.663 


0.00084 


10,644 


0.00084 


10.742 


0.00079 




100 



200 



300 



400 



500 



10- 



Pigure 1.- Velocity distribution in the absence of heat transfer. 
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Pigore S.- felocity distribation in the presence of heat ' 
transfer, for small velocities, m • 
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Figure 2.- Velocity distrilsution in the presence of heat 
transfer, for small velocities. 
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Figure 4a.- Dependence of the drag coefficient c^, referred to 

the density IS in the free stream, on log for 
(JU = 0 (absence of heat transfer but with compressibility effect). 
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Pigure 5.~ Dependence of Cf on lo^ R^^ for 

U) = 0.05 (at small heat transfer). 

(Reference to p in the free stream) 
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Figure 6.- Dependence of c^ on log 11^ for 

cu =0.1 (for large heat transfer), 

(Reference to p in the free stream) 



